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Abstract 



Nonrelativistic QED bound states are difficult to study because of the pres- 
ence of at least three widely different scales: the masses, three-momenta (p, 

oo ; 

and kinetic energies (Ki) of the constituents. Nonrelativistic QED (NRQED), 

On 



an effective field theory developed by Caswell and Lepage, simplifies greatly 
bound state calculations by eliminating the masses as dynamical scales. As 



we demonstrate, NRQED diagrams involving only photons of energy Ey ~ p$ 
contribute, in any calculation, to a unique order in a. This is not the case, 

5_J ■ however, for diagrams involving photons with energies E-y ~ Ki ("retarda- 

d _ 

tion effects"), for which no simple counting counting rules can be given. We 
present a new effective field theory in which the contribution of those ultra-soft 
photons can be isolated order by order in a. This is effectively accomplished 
by performing a multipole expansion of the NRQED vertices. 
It is remarkable that the spectrum of the hydrogen atom is one the first application 
of quantum mechanics being taught and yet it is almost never mentioned in textbooks on 
quantum field theory and QED. Even when the problem of bound states is mentioned, it 
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is made clear that it is a difficult subject and that, to quote Ref. |L|, "accurate predictions 
require some artistic gifts from the practitioner" . 

The problem in studying bound states with relativistic quantum field theory is that the 
conventional perturbative expansion in the number of loops breaks down completely. The 
physical reason is the presence of energy scales absent from scattering theory. Indeed, the size 
of an atom made of two particles of charge — e and Ze is of order the Bohr radius ~ 
(where \i is the reduced mass) which, by the uncertainty principle, provides an additional 
energy scale ~ Z^ia. Because of this new energy scale, there is a region of the momentum 
integration in which the addition of loops will not result in additional factors of a. Moreover, 
if Z <C 137 (condition to which we restrict ourselves in this paper), this energy scale is much 
smaller than the masses of the particles and the system is predominantly nonrelativistic, a 
simplification not taken advantage of in traditional approaches. In addition, a third energy 
scale, again vastly different from the previous two, is set by the particles kinetic energies 
(~ {Z jia) 2 / rrii) and further complicates bound state calculations. 

The problem is greatly simplified by using a Schrodinger theory corrected by the usual 
relativistic corrections obtained by performing a Foldy-Wouthuysen-Tani || transforma- 
tion to the Dirac lagrangian and expanding in powers of p/m. These corrections include 
the Darwin interaction, the spin-orbit coupling, the relativistic corrections to the energy 
(— p 4 / (8m 3 ) + p 6 / (16m 5 ) + . . . and so on (from now on, it is this corrected theory that we 
will refer to as the "Schrodinger theory"). The effects of these interactions can be computed 
by applying Rayleigh- Schrodinger perturbation theory using the Schrodinger wavefunctions 
as unperturbed states, a process familiar from elementary quantum mechanics. Calculations 
are much simpler in this framework both because it takes advantage of the nonrelativistic 
nature of the problem and because it sums up the effects of the Coulomb interaction, re- 
sponsible for the breakdown of covariant perturbation theory, into the wavefunctions. 

However, such a theory is useless for high precision calculations. This is because it does 
not contain the physics corresponding to the high energy (p ~ m) modes of either the 
fermions or the photon. This has two consequences. The first one is the appearance of 
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divergent expressions. These divergences show up in second order of perturbation theory 
(PT) as well as in first order of PT if sufficiently high order (in 1/m) operators are considered; 
for example, the operator p 6 /(16m 5 ) mentioned above is divergent when evaluated in first 
order of PT. These divergences are due to the fact that this theory reproduces faithfully 
QED only when the momenta probed by the interactions are much smaller than the electron 
mass. This condition is not satisfied in most interactions considered beyond first order PT, 
or when the operators contain sufficiently high powers of derivatives to probe the relativistic 
(p ~ m) behavior of the wavefunctions. 

The second consequence is the absence in the Schrodinger theory of operators correspond- 
ing to QED diagrams with photons of energies ~ m e , such as the process e~e + — > 7 — > e~e + 
and the decay of of an electron-positron pair into an odd or even number of photons. These 
processes are clearly important; the first contributes to the lowest order hyperfine splitting 
in positronium, and the others cause the decay of the ortho (total spin S =1) and para 
(S = 0) states of positronium. 

Let us emphasize again that these problems are due to the fact that, in a quantum field 
theory, the high energy modes cannot be simply discarded; they play an important role, 
even in processes involving only nonrelativistic external states. 

Caswell and Lepage ( @, @) have shown how to modify the Schrodinger theory to 
incorporate relativistic effects in a consistent and rigorous manner. They constructed an 
effective field theory (eft) that reproduces QED in the nonrelativistic regime ( p <C m e ) and 
which they christened nonrelativistic QED (NRQED). Although NRQED has been around 
for more than ten years and have been used in high precision calculations in positronium and 
muonium ( 0, [|], ||) , it is still little known, both by the atomic physics community and by 
the eft afficionados. Indeed, the Euler-Heisenberg lagrangian, which describes the scattering 
of photons at energies much below the electron mass, is still the conventional example cited 
as an application of efts in the context of QED. However, the Euler-Heisenberg lagrangian 
(which is a subset of the NRQED lagrangian) has a range of applications quite limited 
which does not include, in contradistinction with NRQED, the important topic of bound 
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state physics. 

In the next section we will review the construction of NRQED. As any eft, NRQED 
contains an infinite number of interactions and is therefore nonrenormalizable. This is not 
a problem because an effective field theory is to be used within a restricted range of energy 
(p -C m e in the case of NRQED) so that only a finite number of interactions will contribute 
to any given process. Which interactions are to be kept for a given precision (in a) is 
dictated by counting rules which are an essential ingredient of any eft. The counting rules 
of NRQED are one of the focus of this paper. 

Clearly, NRQED can be applied to both low energy scattering and nonrelativistic bound 
states. In applications to bound states, the NRQED counting rules are more involved than 
in most eft's because of the presence, as noted above, of more than one dynamical scale in 
the theory: the fermions three-momentum ~ Zfia, and their kinetic energies ~ (Zfia) 2 /mi. 
For the sake of conciseness, from now on we will refer to these two scales as, respectively, 
the "soft" and "ultra-soft" energy scales, E s and E us . Because of the the presence of these 
two scales, there is, in general, no simple connection between an NRQED diagram and the 
order (in a) at which it contributes. 

In this paper, we show how to disentangle the contributions from these two scales in 
such a way that each diagram will contribute to a unique order in a. The first step is well 
known and relies on time ordered (or "old-fashioned") perturbation theory together with 
the Coulomb gauge to separate the "soft" photons (with energy E y ~ E s ) from the "ultra- 
soft" ones (Ej ~ E us ). The counting rules for the diagrams containing only soft photons 
are straightforward and a one-to-one correspondence between a diagram and the order of its 
contribution can be established. The diagrams with ultra-soft photons are more complicated; 
not only do they contribute to an infinite number of contributions of different order (in a) 
but in addition the lowest order is not given in terms of simple rules. This leads us to the 
second step in our separation of scales, which amounts to performing a multipole expansion 
of the vertices involving ultra-soft photons, leading to a new (infinite) set of interactions. 
This can be interpreted as defining a new effective field theory which is superior to NRQED 
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for dealing with ultra-soft photons and which we will call "MQED" (for "multipole QED"). 
We will show that in MQED, diagrams containing ultra-soft photons contribute to a unique 
order in a, as given by new counting rules. In addition to having simple counting rules, 
MQED is better adapted to the study of processes involving ultra-soft photons such as the 
Lamb shift or the generation of certain types of logarithms. These topics will be addressed 
elsewhere @. 

Our paper is divided as follows. In section I we introduce NRQED and its Feynman rules, 
in the context of time ordered perturbation theory. In section II we show how time ordered 
PT permits to separate the contributions from soft and ultra-soft photons, and give the 
counting rules for diagrams containing only soft photons. In section III, we first illustrate 
the breakdown of the previous counting rules for diagrams containing ultra-soft photons. 
We then show how to incorporate the multipole expansion in the NRQED vertices, and how 
this leads to defining a new theory which we will call "MQED" for "Multipole QED". In 
section IV we give the general MQED counting rules and some examples to illustrate their 
use. 
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I. NRQED 



In principle, there are two ways of deriving an effective field theory if the underlying 
theory is know. Firstly, one can integrate out the modes of energies > A p h ys where A p h ys 
is the energy below which the effective theory is to be used (we will keep the subscript 
phys to distinguish this A from the regulator cutoff to be introduced later on; in NRQED, 
A p hys — m). In practice this is technically difficult to do or even impossible, as in the case 
of low energy QCD. The second method consists in writing down the most general effective 
field theory composed of the low energy fields and consistent with the symmetries of the 
underlying theory. The eft is not restricted by renormalizability and contains therefore 
an infinite number of operators, each accompanied by an independent coefficient. If the 
underlying theory is perturbative in the range of energy E < A p h ys , then these coefficients 
can be computed, order by order in the loop expansion, by setting equal, or "matching", 
some scattering process computed in both the underlying and the effective theories. In the 
case of low energy QCD, where such a matching is not possible, the coefficients must be 
determined phenomenologically and the usefulness of the eft is restricted by the wealth of 
data available. 

For NRQED, we follow the second method which requires to first identify the low energy 
degrees of freedom and the relevant symmetries. There will be a field for the photon and 
one for each of the charged particles participating to the process under study such as the 
electron, the positron, the muon, proton, etc. Notice that the fermion fields correspond to 
two-components Pauli spinors. A particle and its associated antiparticle are independent 
fields in a nonrelativistic field theory; they simply correspond to distinct particles of opposite 
charge. NRQED must obey the symmetries of low energy QED such as invariance under 
parity, Galilean and gauge invariance, etc. Lorentz invariance is not necessary except for 
the terms containing photon fields only. 

It is convenient to decompose the NRQED lagrangian in the following way: 

£nRQED = ^2-Fermi + £-4-Fermi + ^-photon + • • • (1) 
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where C 2 -Fermi and C^-Fermi are the interactions containing two and four fermions, , re- 
spectively, and Cphoton is the pure photon lagrangian which includes the Euler-Heisenberg 
lagrangian. We will not display the operators containing six or more fermions fields which, 
in all practical applications, can be ignored because their contribution is suppressed. The 
lagrangian C 2 _ Fer mi is given by 

D 2 D 4 

C-2-Fermi = W { iD t + o h o i + Ci<7 ■ B 

+c 2 (D • E - E • D) + c 3 cr • (D x E - E x D) + . . (2) 

where ip represent a two-component fermion field of charge q and mass m. Notice that in 
NRQED, a particle and its associated antiparticle differ only by their charge. 
The first few terms of C^-Fermi are given by 

c 4 ^(T(-ia 2 )(x f ) T ■ X T (i°2)erip+ c 5 ip\-ia 2 ){x^) T X T {^2)1^ 

+ CQ(^l)\-ia 2 )(jT) 2 x ■ xH^2)crip + h.c.J + c 7 ip^ 'crip ■ X 1 + c sV ; Vx t X + --- (3) 

The first three terms are only present when tp and \ are associated with a particle and its 

antiparticle such as the electron and positron; they come from QED annihilation diagrams 

(the factors of a 2 and the transverse operator T are necessary because we are using the same 

definition for both the particle and antiparticle spinors, see dig)). 

As will become clear in our derivation of the counting rules, the Coulomb gauge is the 

most efficient gauge for the study of nonrelativistic systems. In this gauge, the first few 

terms of C phot0 n are @ 

1 k 4 k 4 k l ki 

- -F, U F^ + c 9 A°(k)—A°(k) - c 9 A\k)— A\k) (5 l3 -_£-) + ... (4) 
4 m z m z k z 

Before discussing the calculation of the coefficients q, we will switch from the lagrangian 
to the hamiltonian. We do so because the counting rules in a nonrelativistic bound state 
are most easily derived in the context of time ordered (or "old-fashioned") perturbation 
theory (TOPT for short) and in TOPT one must work with the hamiltonian rather than the 
lagrangian. We remind the reader that, in contradistinction with covariant PT, in TOPT 
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the vertices conserve only three-momenta and the virtual states are always on-shell. The 
total energy, however, is not conserved by the intermediate state so that, in this formalism, 
it is the violation of energy that characterizes the virtual state rather than the off-shellness 
of the particles, as in covariant PT. 

Using D = z(p — qA) and D t = d t + iqA , the NRQED hamiltonian is given by 



7~{-2 -Fermi — "0 



^ + ^o-/ i -^-(p , + p)-A+^A-A 

2m Sm 6 2m 2m 

- ida • (k x A) - c 2 k 2 A° 

+ 2c 3 <x ■ (p' x p)A° - 2qc 3 a ■ (k a x A(ki))A°(k 2 ) 



+ c 3 k°a ■ ((p' + p) x A) + 



V'(P) + XX terms (5) 



Wi-Fermi = ~c 4 ip ] (r(-ia 2 )(x ] ) T ■ \ T {^2)0^ - c b 4) ] {-ia 2 ){x ] ) T X T {iv2)ip + ■ ■ ■ (6) 

1 k 4 k 4 icty 

n photon = -(E 2 + B 2 ) - c 9 A°(k)—A°(k) + c^i^—A^k) {5 i3 --_) + ... (7) 

A 771 771 rC 

As explained previously, the coefficients are determined by computing some low energy 
scattering process in both QED and NRQED and matching the results. The coefficients 
of the operators in Ti.2-Fermi can be computed by considering the scattering of a charged 
particle off an external field (see Ref. or || for an explicit matching). The coefficient 
C4 is obtained by matching the tree level QED annihilation diagram e + e~ — > 7 — > e + e~ to 
the NRQED interaction. The tree level contribution to C5 comes from the QED diagram 
e + e~ —>■ 77 —>■ e + e~ and is therefore of order a 2 . On the other hand, eg comes from the 
one-loop vacuum polarization. One finds 

Q 



2 



2m 8m 

iq air 

C3 — — C4 — — 

8m 2 m 2 

a 2 a 
c 5 = — (2-21n2 + m) c 9 = — . (8) 
m l 157T 

The imaginary part of C5 corresponds, via the relation Im(E) = — T/2, to the decay rate of 
positronium in a singlet (S = 0) state, the quantum number carried by the corresponding 
operator. 
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Notice that the relation between the powers of a and the number of loops is broken in 
NRQED, since factors of the coupling constant arise from coupling to all photons whereas 
the eft contain only photons with momenta |k| <C m. For the same reason, the tree level 
matching, for example, involves NRQED tree diagrams, but may involve QED loop diagrams. 
By "tree level matching" we will mean matching involving tree level NRQED diagrams. 

The one-loop matching modifies the values of the tree level coefficients so that we will, 
from now on, write the coefficients in the form 

c H — > d5i (9) 

with 5{ = 1 + O(a). As in conventional renormalization, tree level as well as one-loop 
NRQED diagrams enter in the one-loop matching and this defines the 0(a) corrections 
to the NRQED parameters; the only difference with conventional renormalization is that 
the calculation is matched to a QED result instead of an experimental input. Because 
the one- loop NRQED integrals are divergent, they must be regularized. There are many 
possible regulators; one can use dimensional regularization or a simple cutoff A^j on the 
momentum integrations (which is permitted because NRQED breaks Lorentz invariance to 
start with). The NRQED coefficients defined by the matching are then cutoff dependent, 
i.e. they must be viewed as bare parameters. In contradistinction with QED, the divergent 
terms are not only logarithmic but power-law, (A#/m) n , as well. This cutoff dependence 
is of course canceled in any physical calculation, by invariance under the renormalization 
group. Obviously, one can also set A R = A phys = m directly, but since the bare coefficients 
are then finite, this can be misleading if one is not careful about renormalizing the effective 
theory properly (for a more thorough discussion of this point, see ||). 

The one-loop matching of some of the coefficients appearing in (0) has been performed 
in Refs. B and |10[ and the corresponding S^s appearing in (0) were found to be 



6 1 = 5 F = 1 + a e + 0(a 

a 



S 2 = S D = 1 + -- 

7T 3 



. TYl . 11 

'"W + 24 



+ 2a e + 0{a 2 ) 



5 3 = 5s = l + 2a e + C(a 2 ) 
44a 

5, = 5 A = 1 - = (10) 

97T 

where a e is the electron anomalous magnetic moment which, to the order of interest, can be 
taken to be a/(27r). We have redefine our coefficients to follow the convention of (but 
notice that our 5 correspond to their coefficients c); the subscripts F, D, S and 4 — F stand 
for Fermi, Darwin, spin-orbit and four-Fermi interaction, respectively. 

We now turn to the task of writing a general form for the NRQED coefficients. Before 
doing so, we must address the issue of the photon mass, which provides an additional scale 
and has the potential of complicating our analysis. The photon mass does not appear in 



(1C) but this might appear fortuitous. However, since any photon mass dependence is a 
sign of sensitivity to very low momenta and NRQED is designed to be equivalent to QED 
in this region of phase space, any infrared singularity in a QED diagram is also present in 
the corresponding NRQED diagram, so that it gets canceled in the matching. Therefore, in 
general, the NRQED bare coefficients do not depend on the photon mass, to any order in 
the matching. From this, it follows that the coefficients have the general structure 

oo 

c i (A R ,m 1 ,m 2 ) = c° a Ui 5 i {A R ,m 1 ,m 2 ) = c° a n * [1 + £ a 1 ' c^(A K , m u m 2 ) ) (11) 

k=i 

where is now a generic symbol representing any NRQED coefficient and the ij on the 
coefficients q is an index, not an exponent. We have decomposed the lowest order term as 
a coefficient c° of order one times a factor a Ui which is different for different operators. As 
an example, the Darwin interaction, which contains the factor Ci, has n\ = 1/2 whereas 
the singlet annihilation operator, which contains C5, has n.5 = 2. The index U indicates the 
number of loops used in the matching. 

The coefficients c contain, in general, finite pieces plus power-law terms as well as loga- 
rithms divergent terms. Notice that for a fixed /», there are, in principle, an infinite number 
of terms to calculate because there are an infinite number of Zj — loops NRQED Feynman di- 
agrams but only a finite number of interactions must be considered in any given calculation, 
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as specified by the counting rules (which will also dictate the order at which the matching 
must be performed). 

The Feynman rules for the first few interactions of Eqs.(|7|), ((^) and (||) are given in 
Fig.[l]. We will draw the diagrams with the time flowing to the right. In the rules for 
the vertices we have followed the example of || and used the expression "dipole vertex" to 
represent the p • A interaction even though, as pointed out in || , the NRQED hamiltonian 
is not an expansion in multipoles. Also, we have use some Fierz reshuffling to rewrite the 
annihilation vertex in the form given in Fig. [1]. As for the propagators, we have used time 
ordered perturbation theory where there is one propagator for each different intermediate 
state, defined by cutting the diagram with a vertical line. The general rule for a time-ordered 
propagator is 

1 



Eq — Ei 

times a factor 



(12) 



for each transverse photon present in the intermediate state. In (|T2D, E stands for the energy 
of the initial state and Ei for the energy of the intermediate state. One uses nonrelativistic 



energies, p 2 /(2m), for the fermions and Vk 2 + A 2 for the photons. In Fig. [1] the propagator 
is given for an intermediate state containing only one fermion or one transverse photon. In 
Fig. [2] , the corresponding expressions are given for the states containing two fermions or two 
fermions plus one transverse photons, which are the situations most often met in NRQED 
calculations. 

One must sum over all the possible time ordered diagrams and integrate over all the 
internal three-momenta, with a measure d 3 p/(2n) 3 . Notice that we prefer to include the 



factors of l/(2y/k. 2 + A 2 ) corresponding to the transverse photons in the propagators instead 
than in the measure for reasons that will become clearer below. 

In this work we will be mainly interested in applications of NRQED to bound state 
calculations in which case the external lines are not associated with free spinors but with 
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wavef unctions. In general, the wavefunctions are obtained by solving a Bethe-Salpeter type 
equation, with some approximated kernel. This is equivalent to summing up an infinite num- 
ber of this kernel into the wavefunctions. We will show below that the NRQED counting 
rules single out (in the Coulomb gauge) the Coulomb interaction as being the only nonper- 
turbative interaction in a nonrelativistic bound state so that this part of the analysis reduces 
to solving the usual Schrodinger equation. In our explicit examples, we will use the ground 
state wavefunction, given by 



where 7 = Z\ici (the energy of the state is given by — 7 2 / (2/x) ) and £1, £2 are the spinors of 
the two particles making up the bound state, with 



We will not write down the states of higher angular momentum since they are, for the 
purposes of establishing the counting rules, equivalent to the above states (for the momen- 
tum Schrodinger wavefunctions for arbitrary quantum numbers, see ||11||). As just men- 
tioned, using Schrodinger wavefunctions for the external states means that we are summing 
the Coulomb interaction between the external legs. All other interactions can be treated 
perturbatively, which will be shown to be self consistent with the counting rules. 



We now consider a nonrelativistic bound state made, to simplify the discussion, of two 
particles of equal masses and of charges ±e. We will also assume that it is in its ground 
state (n = 1). We will generalize our results at the end of this section. 

There are two important energy scales in such a bound state, the typical bound state 
momentum 7 and the binding energy — 7 2 /m. For a nonrelativistic fermion, for which the 
dispersion relation is given by the usual E = p 2 / (2m) , using either scale leads to pf ermion ~ 7 




(14) 




(15) 



II. COUNTING RULES: SOFT PHOTONS 



12 



(from now on, by p and k we will always mean the magnitude of three-momenta). In the case 
of a photon, for which E = k, using the bound state momentum or binding energy yields two 
very different scales for k, namely k ~ 7 = ma and k ~ r y 2 jm = ma 2 . We will refer to this 
first type of photons as "soft" photons, and to the second type as "ultra-soft" photons. For 
the sake of completeness, we define "hard photons" as the photons with fc~mor greater. 
These photons play no dynamical role in NRQED since they have been integrated out of 
the theory and their only effect is buried in the theory's coefficients. 

The first step in deriving counting rules is to separate diagrams involving soft photons 
from diagrams with ultra-soft photons, since they bring in very different scales, which will 
necessarily complicate the rules. This is where the use of the Coulomb gauge in conjunction 
with time ordered PT will be crucial in simplifying the analysis. 

Consider a transverse photon exchange between two fermions in a nonrelativistic bound 
state. This is represented by the two time ordered diagrams of Fig. [3], where we put the time 
axis toward the right and the \1/ attached to the external lines represent the wavef unctions. 
The photon will contain both soft and ultra-soft components. Now, if the photon is soft, its 
momentum as well as its energy are of order the fermion momentum ma so that its energy 
is much greater than the fermion energies. This means that from the point of view of the 
fermions, the propagation of the soft photons is instantaneous and is therefore represented 
by vertical lines in time-ordered diagrams. 

This can be seen more qualitatively by looking at the explicit expression for the inter- 



mediate state propagator, which is given by (recall that k = y |k 2 | 



2k v tJ k 2 ' V-7 2 /(m) - p 2 /(2m) - (p - k) 2 /(2m) - k 



+ - 7 2 /(m) -p 2 /(2m) - (p-k) 2 /(2m) - k)' ^ 
Notice that the photon mass can be set to zero in bound states calculations, the size of 
the atom preventing the appearance of any infrared singularity; the scale of the fermion 
three-momentum p is of order 7. For soft photons, k ~ Zfia, and we clearly see that k 
dominates in the propagators for the intermediate state so that we can approximate (|16|) by 
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which corresponds to a single diagram, with an energy independent photon propagator. This 
corresponds to the transverse photon propagator of || if one approximates fcjj — k 2 m — k 2 
and set A = (this is why we kept the l/(2k) factor in the definition of the propagator 
instead of the measure). This shows again that in a time ordered diagram the propagation 
of such a photon is represented by a vertical line, i.e. an instantaneous interaction, since it 
is independent of k so that its Fourier transform contains a delta function in time. 

We can now isolate the soft from the ultra-soft components in any photon exchange 
by rewriting the time ordered diagram as a sum over an instantaneous interaction and a 
"retarded" one, as in Fig. [4] (this is why we refer to the effects of ultra-soft photons as 
"retardation effects"). If we restrict ourselves to NRQED diagrams containing soft photons 
only, then all photon exchanges can be represented by vertical lines. In real space, such 
interactions are represented by potentials local in time, i.e. functions of |r x — r 2 | only. 

Besides photon exchanges, the only other possible interactions are the self energy inter- 
actions such as —p 4 /8m 3 and the contact interactions contained in C^-femn, C>6-Fermh etc 
These can clearly be represented by potentials, so that an arbitrary diagram containing soft 
photons only can be written as a string of potentials connected by fermion lines only. In this 
case, the intermediate states contain fermion lines only and the time ordered propagators 
take on a particularly simple form. If there are no interaction from £ n >4, for example, the 
propagators are all of the form 

1 = 1 = m (18) 

Eq — E (intermediate state) — 7 2 / m — p 2 / (2m) — p 2 / (2m) 7 2 + p 2 

where we have used the fact that the bound state energy is — 7 2 /m. The crucial observation 
is that the mass dependence factors out, in the form of the overall factor of m, leaving 7 
as the only dynamical scale in the integrals. If interactions contained in £ n> 4 are included, 
then some intermediate states will contain more than two fermion lines, but it will always 
be of the form 
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1 



m 



(19) 



7 2 /m-£ iPl 7(2m) 



2 + Ei P 2 



and the mass still factors out. Therefore any NRQED diagram containing only soft photons 
leads to an integral of the form 



where b is some integer that depends on the types and number of potentials. The product is 
over all the vertices of type j, with coefficients Cj, as given in (|TTD . Again, the crucial point 
for the following discussion is that the mass m does not appear in the integrand, i.e. does not 
play any dynamical role. There are two scales in the integral, 7 and A#, but the invariance 
under the renormalization group implies that the divergent Ar dependent terms arising from 
the integrations will be canceled by corresponding terms in the bare coefficients c^(A#). As 
noted before, these divergent terms are either power-law, i.e. of the form (A#) n with n being 
a positive integer, or logarithmic. The power law terms are canceled exactly whereas the A^j 
in the logarithms get canceled after combining logarithms containing different scales which 
leaves, in the end, logarithms of a. 

How the logarithms become finite is instructive in that it clearly illustrates the separation 
of scales accomplished by the effective theory. As mentioned in section 2, some NRQED 
bare coefficients contain divergent logarithms of the form ln(A^/m) (as is explicit in (|TTD ). 
To be precise, the NRQED scattering diagrams appearing in the matching process contain 
logarithms of A# over m since these are the only two dynamical scales of the eft, whereas 
the QED scattering diagrams contain logs of the form ln(ra/A); upon solving for the bare 
coefficients, the logarithmic dependence is then of the form ln(A#/m) (again, the photon 
mass dependence drops out entirely for the reasons explained above). On the other hand, the 
NRQED bound state integrals can only depend on the scales A# and 7, yielding ln(Aj?/7). 
In the end, the logs of the bare coefficients combine with the logs generated by the bound 
state integrals to give corrections of the form ln(7/m) = In a. We see how the use of an 
effective theory has separated the contributions from all the scales present in the problem 




(20) 
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(A, 7, m and Ar) in such a way that only two of them played a dynamical role in any 
given stage of the calculation (QED and NRQED scattering diagrams in the matching and 
NRQED bound state diagrams). 

The only A# dependence remaining is therefore of the form (7/A#) n which, upon setting 
Ar = m, leads to corrections beyond the order of interest; in analytical calculations one 
can get rid of these terms by simply letting Ar — > oo at the end of the calculation, as 
in conventional renormalization. In numerical calculations, one needs to keep Ar finite 
because intermediate steps of the calculations are divergent (unless bound state diagrams 
and counterterms are combined before carrying the integals); one must then restrict Ar to 
the range 7 <C Ar Cm in order to assure convergence of the expansion in 1/m used in the 
effective theory. In this case one must keep track of the finite corrections due to the finite 
value of the cutofiQ. 

It is now a trivial matter to write down the counting rules for an arbitrary bound state 
diagram containing only soft photons, i.e. the order in a at which it will contribute. There 
are two sources of factors of a. First, there are the explicit factors contained in the NRQED 
vertices. Secondly, there is a factor of a for each factor of 7 generated by the diagram. To 
be more rigorous, the factors associated to the vertices are genuine factors of the coupling 
constant whereas the factors of 7 are associated with factors of v which scale is set by the 
bound state to be of order a; here it is not important to distinguish between the two types of 



Obviously, dimensional regularization can be used instead of a momentum cutoff. The power law 
divergences are then either entirely absent or replaced by 1/e divergences. Again, these divergences 
cancel, by invariance under the renormalization group. This leaves logarithms depending on the 
scale /i, which gets canceled in the way described above for the logA^ terms. In actual explicit 
analytical calculations, using dimensional regularization or a momentum cutoff is simply a matter 
of taste. However, for high precision calculations, where numerical calculations are required, an 
explicit cutoff is necessary. 
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contributions but this is necessary in QCD bound states because of the noticeable running 
of the strong coupling constant [|12|], ||IBf . 

By simple dimensional analysis, there will be a factor of 7 to compensate each explicit 
factor of mass appearing in the vertices and each factor of mass due to the fermion pair 
time ordered propagators. An arbitrary bound state diagram is built out of a given number 
of potentials, Np, connected by Ntop time ordered propagators. For example, consider 
Fig. [5], where 3 potentials are connected by 2 time-ordered propagators so that Np = 3 
and Ntop — 2 for that diagram. For later use, we also define Vj as the number of vertices 
contained in the i th potential and V as the total number of vertices in the diagram 

Np 

V = 5>i- (21) 

We now define the "vertex mass degree" dj as the number of inverse masses contained in the 
jth ver ^ ex aric [ "potential mass degree" T>i as the number of inverse masses contained in 
the i th potential, 

V t = Y^d r (22) 

For example, the first potential of Fig. [5] (the self energy potential) has T>\ = 3, whereas 
T>2 = (for the Coulomb interaction) and = 4. Since each potential generates T>i factors 
of inverses masses and each fermion-fermion time ordered propagator generates one factor 
of m, an arbitrary diagram having the dimensions of energy will then generate a factor 7 A , 
with 

v 

A = 1-N T op + Y,^ (23) 

i=l 

where the sum is over all the vertices in the diagram. For the present purposes, it is more 
convenient to write A as 

X = l-N TOP + J2^ (24) 
where now the sum is over all the potentials in the diagram. 
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We know define the "coupling constant degree" C, as the total number of explicit factors 
of a contained in each potential, namely 

C t = + ^i) (25) 

j 

where rij and are the powers of a associated with the coefficient of each vertex as defined 
in©- 

Finally, a diagram made of Np potentials will contribute to order ma^ with ( being the 
sum of Eq. fl24|) and the coupling constant degrees ( [25]) of all the potentials: 

N P 

( = J2(V i + C t ) + l-N TOP . (26) 
i=i 

It is easy to see that Ntop an d N-p are related by Ntop = N-p — 1 so that we can write 

N-p N P 

( = J2(V 2 , + C t )+2-N P = £(V i + C i -l) + 2. (27) 

i=l i=l 

Ths expression gives the order in a of any NRQED diagram containing only soft photons, 
keeping in mind that this result can be enhanced by factors of ln(a). For the example of 
Fig. [5], one finds ( = 8. 

Eq.fl27p shows clearly that if there is a potential for which T>i + Cj = 1, perturbation 
theory will break down and it will have to be summed up to infinity. It is an easy matter to 
find such a potential. We can choose lj = (i.e. the coefficients of the vertices have their tree 
level values). Now, rij is zero for the self-energy interactions, but the lowest value that the 
mass degree can take is 3, corresponding to the interaction — p 4 /(8m 3 ), so that the condition 
T>i + Ci = 1 cannot be fulfilled. Many potentials have rij = 1 (i.e. one factor of a) but the 
only one with, in addition, T>i = (no inverse masses) is the Coulomb potential — e 2 /k 2 . 
Therefore, as expected, only the Coulomb interaction must be summed up to infinity and 
the resulting contribution is, from (pTD, of order ma 2 ; all other potentials can be treated in 
perturbation theory. 

In an actual calculation, the counting rules are used in the following way. For a given 
process (hyperfine splitting, decay rate, etc.), one selects all the diagrams with the appropri- 
ate quantum numbers that will, using the counting rules, contribute to the order of interest. 
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The counting rules determine not only the diagrams that must be retained, but also, via the 
lj dependence in (|27|), the number of loops that must be used in the matching of each vertex. 
The matching is then carried for each vertex using scattering diagrams in both QED and 
NRQED. For a given number of loops, there are an infinite number of NRQED scattering 
diagrams, but here the counting rules are used a second time to pick the NRQED diagrams 
that need to be taken into account. Notice that in the matching process, which involves 
scattering diagrams, one uses ([27]) even though this relation was derived for bound state 
diagrams. Once all the relevant diagrams have been taken into account and the NRQED 
coefficients have been renormalized to the appropriate order, the final calculation will be 
finite and will reproduce the QED result, to the order of interest. 



A. Extension to arbitrary masses and charges 

We now extend our counting rules for two constituents having arbitrary masses mi 
and m 2 . The above derivation must then be modified at two points. First, the NRQED 
coefficients given by (|11|) will now contain a dependence on mi and m 2 : 



Ci(A R ,m 1} m 2 ) = c i (m 1 ,m 2 ) a n 



1+ a li c l i(A R ,m 1 ,m 2 ) 



(28) 



No simple general expression can be given for the mass dependence of the coefficients c; it 



arises from QED loop diagrams entering the matching and may involve logarithms of mi/m2, 

etc. The mass dependence of the zeroth order coefficients c° can be taken into account in 

the following way: first, define the vertex mass degrees with respect to each mass, dj(mi) 

and dj(m 2 ) as the number of inverse masses m,\ and m 2 contained in the vertex. For a given 

NRQED bound state diagram, one can then define the following two indices 

v 

K = £di(mi), (29) 

i=l 
V 

p = ]T^(m 2 ). (30) 

i=l 

Obviously, such a diagram will contribute an overall factor l/(m1m 2 ). Since the overall 
result must have the dimensions of energy and since the only energy scale provided by the 
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bound state NRQED diagrams is 7, which contains the reduced mass //, the overall mass 
factor will be 

..K+p+l 

W (31) 

The more general counting rule is therefore 

..K+p+l 

O = ^r« C (32) 

times possible factors of ln/ia and functions of the masses mi, rri2 and /i (which, however, 
arise only if some of the NRQED coefficients have been matched beyond tree level). 

Finally, we consider a bound state with constituents of charges — e and Ze. We first 
include a Z dependence in the NRQED coefficients: 



Q(A R ,mi,m 2 ,) = c° i (m 1 ,m 2 ) Z a ' a n ' 



1+ jr a h c l i(A R ,m h m 2 ,Z) 



(33) 



where Oj will denote the explicit power of Z contained in the zeroth order coefficient of the 
jth ver ^ ex Again, the Z dependence of the arises from the computation of QED loop 
diagrams and we will not write a general expression for this dependence but notice that it 
will necessarily be some power of Z. There is an additional Z dependence which, this time, 
we can take into account: an additional Z dependence comes from each factor of 7 = Zfia 
generated by the NRQED bound states. This number is given by (|23|) : 

v 

\ = l-N TOP + '£d i . (34) 
i=i 

A bound state diagram (with all the NRQED coefficients taking their tree level value) will 
therefore generate a factor Z v with rj given by this last expression plus the Z dependence of 
the tree level NRQED coefficients, as given in (p3|): 

v 

r j = l-N T op + ^2(d i + a i ). (35) 
i=i 

Again, the power of Z is independent of the order in perturbation theory for the Coulomb 
interaction since each Coulomb potential increases the sum over cij by one but increases 
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Ntop by one as well. Our most general counting rule for diagram containing soft photons 
is therefore given by 

..K+P+l 

O = Sr^« c ( 36 ) 

m!{m2 

times possible factors of ln(Z/i«), and dependence on mi, m 2 , // and Z arising from the loop 
corrections to the NRQED coefficients. 

III. COUNTING RULES: ULTRA-SOFT PHOTONS 

The above derivation relied heavily on the fact that the only scale present in the bound 
state diagram was 7. However, if we start considering ultra-soft transverse photons, then 
we have to go back to the general time ordered propagator (see Fig. [3]) 

Pa ( 1 



2k V7 2 /(2a0 + (p - k) 2 /(2m!) + p 2 /(2m 2 ) + k 



+ 7 2 /(2/i) + P 2 /(2m l ) + (p - k) 2 /(2m 2 ) + k) (3?) 
where we have defined the transverse projection operator 

Pa - - ^r- (38) 

In general, such a propagator would contain both the soft and ultra-soft scales so that 
counting rules would be impossible to establish. However, we have already isolated the 
soft contribution in an instantaneous interaction with the photon propagator given by ( |T7| ) . 
Therefore, if the contribution from the soft photon is calculated separately, only the ultra- 
soft scale remains in d37|) . We represent this separation graphically in Fig. [4] where a general 
transverse photon (on the lhs) is represented by a slanted wavy line and, on the rhs, the soft 
photon contribution is represented by a vertical (instantaneous) wavy line and the ultra-soft 
contribution is represented by a slanted, broken, wavy line. To get the ultra-soft propagator, 
we must therefore subtract from the general propagator the expression corresponding to the 
soft photon propagator which we have seen in ([IT]) to be — "P^/k 2 (notice, however, that we 
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would operate this subtraction in a diagram like Fig. [9] where there is no corresponding soft 
photon contribution): 



2k W/(2p) + (p - k) 2 /(2m!) + p 2 /(2m 2 ) + k 

+ 7 2 /(2/i) + P 2 /{2m 1 ) + (p - k) 2 /(2m 2 ) + k ~ k) ' (39) 
This expression now corresponds to the propagator of an ultra-soft photon so the scale of k 
is of order pa 2 . Recalling that the scale of p is ~ pa, we can perform a Taylor expansion in 
k/p ~ a. Applying this to ( j39|) gives 



'2k Vlj 



1 ~ + p • k/mi ff (40) 



T 2 /(2^)+P 2 /(2^) + fc fc (V/ M + P 2 / M + A; 



k 2 /(2m0 + (P ■ k) 2 /m 2 (4i) 



7 2 /(2/i) + p 2 /(2/i) + fc) (l 2 /(2p) + p 2 /(2/i) + fc 
— k 2 p-k/mf (k-p) 3 /mf+ 



7 2 /(2/i) + p 2 /(2^) + k ) 3 f 7 2 /(2^) + P 2 /(2/.) + ^ " 



+ 



+ same with mi <-> m 2 (42) 

where the first line contain the zeroth order term plus the first order one (the p • k term), 
the second line contains the second order contribution and so on. 

Since the expansion is in k/p, we expect that each power of k appearing in the numerator 
will be associated with a power of a with respect to the zeroth order term of the Taylor 
expansion (the first term in (|40D). We will show this explicitly for a few terms. 

Consider first the zeroth order propagator. It contains two inverse powers of k, which 
scales like fia 2 , so that it contributes to the counting rules by a factor l/(m 2 a 4 ) (we won't 
distinguish between mi,m 2 and p, to discuss the counting rules). Of course, in an actual 
diagram, other factors will enter to make the overall a contribution of the diagram positive; 
here we are just interested in the relative contribution of the terms in the Taylor expansion. 

Now consider the first order correction (the second term of (fiO"D ) . The numerator k ■ p/m 
scales like ma 2 x ma/m = ma 3 and the denominator contains 3 factors of k so it scales like 
(ma 2 ) 3 . Therefore, the first order propagator scales like 
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ma i /{m i a b ) = l/(mV) (43) 

which is one power of a times the zeroth order propagator. The k 2 /m term in the second 
order Taylor propagator (f|l|) scales like 

(ma 2 ) 2 /(m x m 3 a 6 ) = I /{ma) 2 (44) 

which is down by two powers of a with respect to the zeroth order propagator. It is a simple 



matter to verify that the other term of (f|l]) also contribute with a factor of a 2 with respect 
to the lowest order contribution, and the terms of ( |42f) contribute with a factor a 3 , etc. 

In an actual diagram, the Taylor expansion must of course be carried on the whole 
diagram. As an illustration, we expand the complete integrand corresponding to Fig. [3(a)], 
sandwiched between ground state wavef unctions: 



(2 Pi - h) (kj - 2 Pj ) 



( p 2 + 7 2)2 4 mim2 ^ " W n> 2k 



1 l\ 8v^T 



-7 2 /(2/i) - (p - k) 2 /(2mi) - p 2 /(2m 2 ) - k k) ((p - k) 2 + 7 2 ) 2 



[p 2 + 7 2 ) 4 4 mi m 2 2k (- 7 2 /(2/i) - p 2 /(2^) - k) 2 
(8v^) 2 qi q 2 Vij/ 1 



(p 2 + 7 2 ) 4 4mim 2 2k V (-7 2 /(2/i) - p 2 /(2/x) - Jfe) fc 

1 6gjgj p ■ k 



x 



ApiPj + 2kipj + 2pikj j-^ 3 + . 



(45) 



Again, one can easily verify that each power of k in the numerator is associated with an 
extra factor of a. 

Notice that the spin-spin diagram with an ultra-soft photon, Fig. [8(b)], contains at least 
two powers of k since the NRQED Feynman rule for the Fermi vertex is proportional to k; 
in other words, the first non-vanishing contribution comes from the second order term in 
the Taylor expansion. Therefore, the lowest order contribution of the ultra-soft spin-spin 
exchange is suppressed by two powers of a with respect to the corresponding dipole-dipole 
exchange (this is due to the fact that the Fermi interaction involves the B field). This is 
very different from the corresponding soft photon diagrams which both contribute to the 
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same order. The difference, again, is that only the factors of e and 1/m enter in the soft 
photon counting rules whereas factors of the photon momentum k matter in the ultra-soft 
counting rules. 

Clearly, the fact that one power of a is generated by each term in the Taylor expansion 
will prove crucial in writing down the counting rules of this new, Taylor expanded, version 
of NRQED. However, before doing so, we now want to show that the Taylor expansion we 
just carried is equivalent to a multipole expansion of the NRQED vertices. 



A. Connection with the multipole expansion 

As an example, consider the — q^{p- A + A ■ p) / (2m)ifj interaction contained in the term 
— ip jf 'D 2 / (2m)ip in the hamiltonian. To obtain the NRQED Feynman rule we first replace 
the fields by plane waves: 

p'ry . € ( g-ikTgip-r \ g-ip' r / e -ik r £ . y e «p 



e 

+ qi( + -) (46) 

y V 2m 2m 1 K ' 

where p, p' are respectively the three momenta of the fermion line before and after the 
interaction, and k is the photon three-momentum; e is the photon polarization. Applying 
the derivatives, we get 



Q 
2m 



(k l -2p l >)e^ p ' +k - p > r . (47) 



The exponential leads, as usual, to the conservation of three- momentum p' = p — k (here we 
considered a photon being emitted). Using this to write — 2p = — p — p' — k and discarding 
all factors associated with the external fields, we obtain the Feynman rule 

-fig. w 

The rule is obviously unchanged if we consider an absorbed photon. Now we consider a 
multipole expansion of this vertex, i. e. we expand the photon field 

e ~ ikr = 1 -ik-r + -(-ik-r) 2 + ... (49) 
2 
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In the following, we will use the notation e _lk r = zeroth order + first order + ... to label 
the terms in the multipole expansion. As usual, this expansion makes sense only if kr <C 1. 
The size of r is set by the bound state to be of order the Bohr radius r ~ I/7. For ultra-soft 
photons, we have k ~ j 2 / /a so that kr ~ a and the multipole expansion is valid. Of course, 
it would be nonsensical to use it for soft photons. Also, the multipole expansion is clearly 
the same as the Taylor expansion performed above since the scale of r ~ 1/p. 

We can easily find the rule for the new vertex. Using the first term of the multipole 
expansion, e" lk r = 1 (corresponding to an El transition) in fl46|) , we obtain 

- ^ e - l(p '- p) r (50) 
m 

where now the exponential leads to the condition p' = p, i.e. three momentum is not 
conserved at the vertices when the multipole expansion is used. This can, however, still be 
used to write the rule for the vertex as before, i.e. 

<«> 

Even though the rule for the vertex is the same as before, the condition p' = p simplifies 
greatly the expression for diagrams containing ultra-soft photons and, in particular, the 
propagator. To see this, we first go back to the time ordered photon-fermion pair propa- 
gator ([39|). We now use in addition the fact that the fermion momenta at the vertices are 
unchanged by the emission or absorption of the photon to write ([39]) as 



2k V7 2 / (2/i) + p 2 / (2/i) + k 7 2 /(/i) + p 2 /(2/i) + k k 



Va ( 1 1 



k V7 2 / (2/i) + p 2 / (2/i) + k k 



(52) 



instead of the form (^9|) which was obtained by using p' = p — k. In (|52|), the scale of k is 
set either by 7 2 /(2/i) ~ fia 2 or p 2 /(2/i), but since p is a fermion three-momentum it is of 
order 7, we get in either case k ~ /1a 2 . This shows explicitly that the multipole expansion 
has permitted us to isolate the ultra-soft scale. 
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To obtain the higher order terms in the multipole expansion, one provides a factor 
(±k • Vp)"*/ri! for each vertex connected to an ultra-soft photon, where n is the order 
of interest in the multipole expansion, and a plus (minus) sign is used if the photon is 
absorbed (emitted). In this expression, the gradient must be taken with respect to the 
three-momentum of the fermion line on the right of the vertex. To apply these rules, it is 
therefore necessary to distinguish between the momentum of the fermion before and after 
the interaction, even though we have to set them equal in the end. 

To illustrate this, we will evaluate the first few multipole corrections to the ultra-soft 
photon propagator. Since, as noted above, one must distinguishes the momenta of each 
fermion and the momenta before and after the interaction, we will use the momenta as 
labeled in Fig. [6], with the understanding that one must set 

Pi = Pi = -P2 = -p' 2 = P (53) 

after carrying out the derivatives. Taking this into account, the intermediate state propa- 
gator in Fig. [6(a)] is 

~ 2k V ^k + 72 /(2/i) + p' 1 2 /(2m 1 ) + p|/(2m 2 ) " ^ (54) 
and the propagator of Fig. [6(b)] is 

~ 2k Vij( -k + 72 /(2/x) + p?/(2m 1 ) + p 2 2 /(2m 2 ) ~ k^ (55) 
If we consider Fig. [6(a)], then we only have to consider the multipole expansion of the 
vertex on the upper line since the other vertex will not act on the intermediate state prop- 
agator. We therefore apply, as we did above, the operator — k ■ V p ' on fl54|) to obtain 

Vij -k ■ p/mi 

2A;(fc + 72 /(2 / i) + p7(2 / x)) 2 ' 

In the case of Fig. [6(b)], we apply the operator ik ■ V p ^ on Q55D with, for result, (recall that 

we replace p 2 by — p after differentiating) 

k ■ p/m 2 



2k {k + 7 2 / (2/i) + p 2 / (2fi)Y 
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(57) 



As expected, this is the same as (|5q) with m 2 replaced by m\. The sum of and (|57|) is 
the result of the first order term of the multipole expansion. To be more precise, this is the 
result obtained from considering the first order term in the multipole expansion of either 
vertex. 

The result of the second order multipole can easily be calculated in a similar way. We 
apply the operator (-k • V p ' i ) 2 /2 to (HD and (k • V p ^) 2 /2 to (H) to obtain 

Vuf k 2 /(2m!) (k-p) 2 /m 2 , \ . , 

~ ' 1 1 o + same with m 2 ^m 1 ). 58 

2k \(k + ^/(2fj) + P 2 /(2 M )) 2 (fc + 7 7(2/i)+P 2 /M) 3 V 1 J 

These expressions correspond to keeping the n = 2 multipole term on either of the vertices 
plus the first order term on both vertex, all of which contribute to the same order in a, as 
we will discuss in the next section. 
We also give the third order result: 

Vij ( k 2 k • p/m 2 (k • p) 3 /m? 



. ., — j + same with m 2 — > mi ] . (59) 

2H(Hf/(2/i) + P 2 /(2M)) 3 (A; + 7 V(2/i) + p 2 /(2 /U )) 4 V V ' 

We have recovered the expressions obtained from the Taylor expansion, Eqs. (HO, HT] and 



12). This is not surprising since the Taylor expansion of a function f(x + a) around x = 



can be written as 

f(x + a) x ^ = e x ^f(a) (60) 



d 



and this is what the multipole expansion accomplishes. 

In an actual calculation, the multipole expansion must of course be carried on the whole 
diagram. This is slightly more complex because the wavefunctions must also be written in 
a way that distinguishes the momenta on each fermion line. To illustrate this, we consider 
again the bound state diagram corresponding to Fig. [6(a)] and work out the expression in 
first order of the multipole expansion. We again use the ground state wavefunction ( JH| ) for 
the external states. Taking this into account, the integrand corresponding to Fig. [6(a)] is 
given by 
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%Vkt> qxq 2 , / W / \ Vij 

fi 2 (pj/ mi + p|/m 2 + 7 2 //i) 2 4m im2 lPl + Pl )l lP2 + P2>J 2k 

1 1 8^/7^ : f s 

( -7 2 /(2^) - P?/(2 mi ) - pi/(2m 2 ) - fc + fcV(p'i7™i + V2 2 /m 2 + 7 2 / V) 2 ^ 

The contribution of the zeroth order in the multipole expansion is obtained by simply using 
the relations (|53|) directly in (^l|). The contribution of the first order multipole expansion 
is then obtained by applying on this expression the operator — k ■ V p ' , which is associated 
with the vertex on the left in Fig. [6(a)] plus the operator k ■ V p ^ for the second vertex, and 
then reexpressing the vectors in terms of p using ([53]). The result is 

qiq 2 8TS7rf5) 2 Vij ( 4p;Pjk ■ p/mi 



4mim 2 (7 2 + p 2 ) 4 2k V(-7 2 /(2/i) - p 2 /(2/i) - k) 2 

qiq2 Sy/^rO 2 / 1 1 

A / O ON/1 aJ I O //r-» \ O I / r\ \ T 



4 mi m 2 (7 2 + p 2 ) 4 2A; V-7 2 /(2/i) - p 2 /(2/i) - A; fc 

P 2 + 7 



x ( 2k iPj + 2k iPi - 16 PlPj ^ 2 - P ) . (62) 



This is, as expected, equal to the expression obtained from the first order Taylor expansion 
(fi5|). A similar calculation for Fig[6(b)] gives the same result as ( |5"2"D with mi <-> m 2 . 

Notice that the zeroth order term in the multipole expansion is obtained by setting 
p' = p in the NQRED vertices. In the case of the Fermi vertex, this gives zero since the 
NRQED Feynman rule is proportional to p' — p = k. This means that the first nonzero 
contribution is of the first order in the multipole expansion. Higher order terms are obtained 
as above i.e. by considering the corresponding factor of (±k • V p ) n /n!. 

Even though we have simply recovered the expressions obtained by performing a simple 
Taylor expansion, there is one important reward for doing so: one can use directly the 
Wigner-Eckart theorem and the familiar selection rules derived in quantum mechanics for 
each interaction generated by the Taylor expansion. This has consequences in decays of 
positronium, and in nonrelativistic QCD bound states [T2|]. This will be explored elsewhere 



To summarize, we have seen that, starting from NRQED, separating the soft and ultra- 
soft scales and applying a multipole expansion (or Taylor expanding) the vertices connected 
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to ultra-soft photons generates a new effective theory with its own set of Feynman rules. This 
theory, which we will call "MQED" (for "Multipole QED" ) has the advantage of generating 
bound state diagrams that contribute to a unique order in a. In the last section, we will 
derive the MQED Feynman rules and show some applications of the counting rules. 

IV. MQED COUNTING RULES 

We can now easily extend the counting rules to include diagrams containing ultra-soft 
photons. The concept of potentials is not well-defined, however, when ultra-soft photons 
are present, so we first rewrite the soft counting rule ( J2~3| ) as a sum over vertices instead of 
a sum over potentials: 

N v 

C(soft photons) = V](c?j + rij + lj) + 1 — Ntop (63) 

3=1 

where V is the total number of vertices contained in the diagram. For a diagram containing 
ultra-soft photons, this rule must be changed to 

N v 

C = J^(dj + n 3 + lj) + 1 - N TOP + 2iV 7 + Muitra-soft (64) 
i=i 

where iV 7 is the number of ultra-soft photons in the diagram. The last term, M. u itra-soft can 
be expressed in two different ways, depending on whether one uses a Taylor expansion of 
the diagram or a multipole expansion of the vertices. In the first case, Ai u itra-soft is simply 
the power of k appearing in the numerator. In the second case, M u itra-soft can be written 

as 

M u lt r a-soft = (65) 
i 

where the sum is over the vertices connected to ultra-soft photon and M.i is the order in the 
multipole expansion to which the i th vertex has been expanded. 

Eq.flBip gives the order, in powers of a, at which an arbitrary MQED diagram will 
contribute. The dependence on arbitrary masses is unchanged by the presence of ultra-soft 
photons and is therefore still given by (p2|) . The charge dependence, however, is different 
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when there are ultra-soft photons because the ultra-soft scale is ~ 7 2 //i — fiZ 2 a 2 so that 
the Z dependence is different then in the soft scale 7 = Z\xa. The expression for the charge 
dependence must then be changed from (pop to 



v 

V = 1 - N TOP + Y,( d i + a i) + 2N i + E Mi (66) 

i=l u—sVi 

where, again, the second sum is over the vertices connected to ultra-soft photons only. 
Our final result is therefore that an arbitrary MQED diagram will contribute to order 

. K+p+l 

O = ^-^Z^of (67) 
with ( defined in fl64[) , 77 defined in (|66|) and k and p defined in (|30l) . 



We now give a few examples of the use of (|B~7D. As a first example consider the interaction 
Fig. [7] in hydrogen, where the ultra-soft photon is connected to an electron line. In this 
diagram, di = d 2 = 1 (there is one factor of 1/m on each vertex), % = n 2 — 1/2 (a factor 
e on each vertex), l\ — I2 — (the p • A interaction does not get renormalized) , Ntop — 1j 
A 7 = 1 and, if the zeroth order in the multipole expansion (or in the Taylor expansion) 
is used, M.i = M. 2 = 0. This leads to a contribution of order a 5 . The mass dependence 
is found to be p 3 /ml and the Z dependence is, from (j66|) , Z 4 . This diagram therefore 
contributes to order 



m 2 



a 5 . (68) 



In fact, this result is enhanced by a logarithm ln(Za) and contributes to the Lamb shift. 

Consider now Fig. [8(a)] in positronium so that Z = 1 and m x = m 2 = m e . In this 
diagram, the transverse photon is soft (it is represented by a vertical line). We can therefore 



use (|63|). One has n\ = n 2 = 1/2 and d\ = d% = \. If the tree level expressions are used 
for the coefficients, then this diagram contributes to order m e a 4 . The same diagram will 
contribute to higher order in a if the loop corrections to the coefficients of the Fermi vertices 
are considered (the one-loop correction being, from (|T0|), a/27r). 

As a final example, consider Fig. [8(b)]. Here the photon is ultra-soft. As mentioned 
previously, the first nonvanishing contribution from this diagram contains two factors of k 
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(one from each spin vertex) so that M u itra-soft m (|>|) is at least equal to two. Ntop = 1, 
N 7 — 1 and the other coefficients are as in Fig. [8(a)], if the tree level coefficients are used. 
One the finds that this diagram will contribute to order m e a 7 . 
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FIG. 1. NRQED Feynman rules 
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Feynman rules (continued) 
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Feynman rules (end) 
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Pi 

-> mi Fermion-Fermion Propagator 
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FIG. 2. Time-ordered propagators for two fermions or two fermions plus one transverse photon. 
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p — k 




p 


p — k 










(a) (b) 

FIG. 3. The two time-ordered diagrams corresponding to the exchange of a transverse photon 
(the vertical lines indicate the intermediate states used for the time-ordered propagators). 




FIG. 4. Separation of a transverse photon into a soft, instantaneous contribution (represented 
by a vertical line) and an ultra-soft propagator (represented by the broken wavy line). 
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FIG. 5. Generic bound state potential; the dependence on the charges and on the masses of 
each vertex is indicated. 
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FIG. 6. The two time-ordered diagrams corresponding to a transverse photon exchange with 
the routing necessary to apply the multipole expansion. 
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FIG. 7. Self-energy diagram with an ultra-soft photon. 
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FIG. 8. Spin-spin exchange with a soft photon (vertical line) and an ultra-soft photon. 



U'? I 

— 



FIG. 9. Ultra-soft photon spanning a Coulomb interaction. In such a diagram, one does not 
subtract the soft photon propagator from the intermediate state propagator because there is no 
corresponding soft photon diagram. 
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